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Abstract. 

We present a detailed analysis of the soft X-ray spectrum of the Seyfert 1 galaxy Markarian 509 taken with the 
XMM-Newton Reflection Grating Spectrometer. An underlying power-law continuum and three warm absorber 
phases provide a good fit to the data along with a number of broad and narrow emission lines. Our three warm 
absorber phases each have different ionization parameters and column densities. We identify a low ionization, 
log £=0.89, high outflow velocity phase producing an Fe M-shell UTA feature along with absorption from OVI 
and NVI. There is an intermediate phase, log £=2. f 4, showing absorption from H-like carbon and nitrogen and 
He- like neon and oxygen. The third, high ionization, log £=3.26, low outflow velocity phase contains absorption 
from OVIII, NeX and highly ionized iron. All phases are blueshifted with respect to the systemic velocity with 
flow velocities ranging from -GOkms" 1 to -510 kms -1 . The observed broad emission features have an RMS 
velocity of 8000L3000 kms -1 for the C VI and O VII lines. 
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1. Introduction 

Soft X-ray absorption by ionized ga s was first iden tified 
in Active Galactic Nuclei ( AGN) by iHalpernl l)l984l) and 
recent observations indicate that over half of Seyfert 1 
galaxies may have these 'warm absorbers' (Reynolds 1997, 
George et al. 1998). Since the launch of the XMM-Newton 
and Chandra high resolution grating spectrometers our 
ability to probe soft X-ray spectra of AGN has been 
greatly increased, enabling astronomers to do in-depth 
studies of the X-ray absorbing and emitting photoionized 
gas and to examine how they relate to the system as a 
whole (e.g. Blustin et al. 2005). Gas properties are ex- 
amined to improve our understanding of the kinematics, 
ionization and composition of these systems. The warm 
absorber is often characterised by its ionization parame- 
ter, £. This determines the ionization state of the material 
and is defined as 



e-4 



nr* 



(1) 



where L is the source luminosity (in ergs -1 ), n is the gas 
density (in cm -3 ) and r is the distance of the absorber 
from the centre of the system (in cm) (Tarter et al. 1969). 
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Markarian 509 was classified as a Seyfert 1 galaxy 
bv lKopvlov et al. but is more like a Seyfert 1 / quasi- 

stellar object (QSO) transition galaxy due to its compact 
appearance and high luminosity of 1.3 — 2.6 x 10 44 ergs -1 
in the 2-10 keV range (Weaver et al. 2001). Its close prox- 
imity, z=0. 034397 (Fisher et al. 1995), compared to other 
AGN of similar luminosity, and the low Galactic column 
density in the line of sight (N# = 4.4 x 10 20 cm -2 ; Murphy 
et al. 1996), make it an ideal candidate for studying the 
inner regions of such objects. 

Markarian 509 has well-studied intrinsic UV absorp- 
tion (Kriss et al. 2000; Kraemer et al. 2003) but much 
is still to be learnt from its high resolution X-ray spec- 
tr um. Markaria n 509 w as first identified as an X-ray source 
bv lCooke et al.l l|l978j) and has since been observed with a 
number of instruments including EXOSAT (Morini et al. 
1987), GINGA (Singh et al. 1990), and BeppoSAX (Perola 
et al. 2000). 

IPounds et all l|200lj) studied the XMM-Newton 
Reflection Grating Spectrometer (RGS) spectrum of 
Markarian 509 taken in October 2000. They modelled it 
with a blackbody continuum as well as a weak Fe M-shcll 
UTA region, Ne IX absorp tion and O VII narrow emission 
lines. lYaaoob et ail §003) performed a thorough analysis 
of the Chandra High-Energy Transmission Grating data. 
They modelled the spectrum with two absorption edges 
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from O VII and O VIII along with Ne IX, Ne X and O VIII 
line absorption. They fitted these spectral features using 
one phase of absorption through a photoionized gas with 
a column density 2.06^0 45 x 10 21 cm~ 2 and an ionization 
parameter of log£ = 1.76+ ; 14 . The gas was found to 
be outflowing at 200±170kms _1 and a curve of growth 
analysis determined an approximate turbulent velocity of 
lOOk ms" 1 . 

iKriss et alJ l|2000|) analysed the UV spectrum 
of Markarian 509 obtained with the Far Ultraviolet 
Spectroscopic Explorer. They modelled seven components 
for the UV absorber, which they grouped into two veloc- 
ity phases with outflow velocities kms' 1 and 370 km s -1 . 
They identified one component as having an ionization 
parameter and column density that could be compati- 
ble with an X-ray absorbing gas. With the exception of 
this component, however, they predicted that the UV and 
X-ray spectra are produced in separate zones within the 
AGN. iKraemer et alJ l)2003|) analysed the UV spectrum 
of Markarian 509 using data from the Space Telescope 
Imaging Spectrograph echclle gratings on the Hubble 
Space Telescope. They identified eight components in the 
UV absorber, with a similar range of velocities to Kriss 
et al. They did not find evidence for a component with 
an ionization and column density compatible with the X- 
ray absorbing phase. However, Kraemer et al. concluded 
from the kinematic association that the UV and X-ray ab- 
sorbers may co-exist spatially, for example in the form of 
a two-phase clumpy medium. 

In this paper we investigate the combined RGS spec- 
tral data from two observations by XMM-Newton. In 
Section we describe the observations and data reduc- 
tion, and the spectral fitting is reported in Section We 
discuss the results in Section^and present our conclusions 
in Section |3J 

2. Observations and Data Reduction 

Our data were obtained with the RGS (den Herder et 
al. 2001) onboard the XMM-Newton observatory. This 
spectrometer provides unparalleled sensitivity in the 6- 
38A (0.3-2.5keV) range with an approximately con- 
stant resolution of 70 mA FWHM. The data analysed here 
are a combination of two observations of Markarian 509 
by XMM-Newton. The first took place on 25 October 

20 00 with an exposu re of 23.4 ks and was presented 
bv lPounds et al.1 l(200ll) . The second took place on 20 April 

2001 with an expos ure time of 3 0.6 ks and t he EPIC data 
were presented in IPaee et all l)2003[) and iDadina et alJ 
ll2005h . 

The data were processed using the XMM-Newton 
Science Analysis System (SAS) version 5.2. First and sec- 
ond order spectra from both RGS, and both observations, 
were extracted. To correct for residual artefacts in the 
effective area calibration (Stuhlingcr et al. 2005), the ef- 
fective area of each response matrix was divided by the 
ratio of a powerlaw + Galactic column fit to the XMM- 
Newton rev. 0084 RGS spectrum of the continuum source 



Markarian 421. The Markarian 509 spectra and response 
matrices were resampled and coadded to produce a single 
spectrum and a single respo nse matrix. This method is 
explained in greater detail in IPatre et alJ l|2003t) . the only 
difference being that we do not treat the background sep- 
arately as we use background subtracted RGS files. All 
of the nominal channel boundaries and energy ranges in 
the re-distribution matrix have been multiplied by 1+z in 
order to correct for redshift. To improve signal to noise 
the data were grouped by a factor of 3, resulting in a 
spectrum with 1000 channels ~30mA wide, well sampled 
with respect to the RGS resolution of ~ 70 mA FWHM. 
To correct for Galactic photoelectric absorption, the effec- 
tive area elements of the response matrix were multiplied 
by the transmission as a function of energy for a column 
density of 4.1 x 10 20 cm~ 2 . Finally, the wavelength scale 
was transformed to the rest frame of Markarian 509. In 
this way the data have been corrected for both astronom- 
ical redshift and absorption within our own Galaxy and 
so, from this point onwards, all data refer to the rest frame 
of Markarian 509. The main software used to analyse the 
data was SPEX 2.00, a spectral fitting package created by 
J. S. Kaastra with high resolution X-ray spectra in mind. 

3. Spectral Fitting 

3.1. Initial Analysis of the Data 

Certain regions of the spectrum of Markarian 509 (Fig. ^| 
are extremely complex, such as those between 12A and 
17A and between 20A and 24A. In these areas it can be 
difficult to distinguish absorption from emission features, 
with the risk that some features may be identified incor- 
rectly or missed completely. Therefore, in order to get 
an estimate of the statistical significances of the major 
absorption or emission lines, the data were put through 
a program that identifies all of the significant spectral 
lines present (Page et al. 2003a). Firstly a smooth con- 
tinuum model is determined by averaging with a sliding 
cell, ~ lA wide, and excluding all points more than 2a de- 
viant from the smoothed continuum model (Fig. It 
is seen that the 2a error bound has a number of sharp 
spikes in it, which are due to chip gaps in the instrument. 
Fortunately none of these occur at a point in the spectrum 
where there is an important spectral feature. A grid search 
is then performed adding a Gaussian of varying width, in 
steps of 50kms _1 up to a = 1000 km s -1 , to the contin- 
uum model at each point in the grid and calculating the 
goodness of fit between the model and the data by min- 
imising x 2 (Lampton et al. 1976). This method is justified 
as the average number of counts in each bin is 140. This 
provides a measure of significance for all of the spectral 
features, shown in a graph of A\ 2 against wavelength in 
Fig [2 

The smoothed spectrum gives us a good overall view of 
the continuum shape; from this we are able to see the main 
features present, such as a number of possible broad emis- 
sion lines, the most prominent being centred on ~21.7A. 
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Wavelength (A) 



Fig. 1. The RGS spectrum of Markarian 509 (grey), the 
smoothed continuum used by the line finding program 
(black) and the 2<7 error bounds (dotted), in the rest frame 
of the galaxy. 



We also see a rise in the continuum at ~ 17A possibly 
due to either the O VII edge or an Fe M-shell UTA. We 
now look to the A\ 2 graph where a number of strongly 
significant lines emerge. A A% 2 value of 9 corresponds to 
a 3cr significance for 1 interesting parameter. However, if 
one is scanning in wavelength space and the equivalent 
width of a feature is fitted, then one has to take into ac- 
count the number of trials when assessing the significance 
of any features. Since we have approximately 450 indepen- 
dent resolution elements, we would expect 1.35 statistical 
fluctuations deviating at or more than 3cr. There are three 
absorption lines with greater than a 3tr significance at ap- 
proximately 13. 5A, 19. OA and 33. 7A. These correspond 
to the absorption lines of Ne IX, O VIII and C VI respec- 
tively. There are also some weaker features at 16. 5A, from 
a possible Fe M-shell UTA, and at 24.7A, due to NVII. 
The Galactic O VII absorption line at 2 1.6 A, seen here in 
the Markarian 509 rest frame at 20. 9A, will be removed 
when fitting the spectrum at a later stage. Three emission 
features are also present at 18. 8A, 21. 7A and 22. lA, which 
are likely to be due to O VIII and O VII emission. 

By taking a close look at the wavelengths of all of these 
lines we can calculate a flow velocity by measuring their 
shifts relative to the systemic velocity. These velocities, 
with the 2(7 error bounds, are shown in Table Q along 
with the assumed species and the observed and emitted 
wavelengths. 

3.2. The Continuum and Broad Emission Lines 

Now that we have the basic overall spectral profile and 
some approximate velocities we can go on to study the 
spectrum in more detail. Before looking at the absorption 
profile we first need to get a good fit to the underlying 
continuum. We start by fitting the data with a power law; 
however, there are regions where this is inadequate (for 
X 2 values see Table EJ). A number of broad bumps are 
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Fig. 2. A graph of A\ 2 against wavelength to show the 
significant lines in the spectrum, in the rest frame of the 
galaxy. 

Table 1. Observed (in the rest frame of Markarian 509) 
and restframe wavelengths, and flow velocities derived 
from measuring the wavelengths of the individual lines. 



Species 


A (A) 


A e (A) 


Vel (kms" 1 ) 


NelX absorption 


13.47tg:g 


13.45 


+570l 6 i°o 


OVIII absorption 




18.97 




NVII absorption 


24 74+ 11 

'^-0.04 


24.78 


— 440+ 850 
^* u -240 


C VI absorption 


oo -,+0.02 
OO. 1 l_o.02 


33.73 


-2101^ 


O VIII emission 


18 89+0-02 


18.97 


-2400tJ«j 


OVII(r) emission 


91 79 +0.07 
' ^-0.05 


21.80 


-nootES 


OVII(f) emission 


22 09+ 007 


22.10 


-170l!f 



apparent in the spectrum, as seen in Fig. ^ and we find 
an improved match to the data by fitting a number of 
broad Gaussian emission lines. The spectrum in Fig. ^ 
shows broad peaks at 19. 9A, 21. 7A and 25. 2A. There also 
appears to be a narrow emission feature either side of the 
C VI absorption line at 33. 7A (Fig.[2J), which can be fitted 
with another broad Gaussian line set to the wavelength of 
the C VI line. 

The broad emission lines at 21. 7 A, due to O VII, and 
33. 7A, due to CVI, give an average intrinsic RMS ve- 
locity of 8000± 3000 km s" 1 (Table EJ. We identify the 
25. 2A broad feature as a possible CVI radiative recombi- 
nation continuum (RRC; when electrons recombine with 
the ions, emitting photons). However, the identifications of 
this and the 19. 9A feature are not sufficiently reliable to be 
retained in the spectral model, which comprises a power 
law and two broad Gaussians at 21. 7A and 33. 7A. The 
fitted continuum level at 20A and 25A will remain lower 
than the data, as can be seen from the residual graph in 
Fig. El however, as we shall see in Section 13.41 there are 
no significant absorption features in these regions, so this 
should not affect our analysis of the absorption profile sig- 
nificantly. 
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Table 2. The variation in \ 2 for the different spectral models used; PL = Power law. 



Model 


x 2 


Degrees of Freedom 


Probability 


PL Photon Index 


PL only 






1553 


992 


1.84 x 10" 27 


2 1Q+ 001 


PL + 2 broad Gaussians 






1463 


988 


4.72 x 10" 21 


9 17 +0.01 
' -0.01 


PL + 2 broad + 2 narrow Gaussians 






1416 


984 


3.93 x 10" 18 


9 -,7+0.01 
' -0.01 


PL + 1 XABS components + 2 broad + 2 


narrow 


Gaussians 


1116 


980 


1.52 x 10" 3 


9 1 n+0.01 
z - lo -0.02 


PL + 2 XABS components + 2 broad + 2 


narrow 


Gaussians 


1074 


976 


1.52 x 10 -2 


9 19+O.OI 

^•^-0.01 


PL + 3 XABS components + 2 broad + 2 


narrow 


Gaussians 


1018 


972 


1.51 x 10" 1 


2 14+ - 02 
z - ly: -o.oi 



Table 3. Emission line parameters for both the broad (top) and narrow (bottom) fitted lines. No flow velocity is 
given for the O VII broad emission line because it is a combination of the three lines of the triplet and so an accurate 
blucshift cannot be calculated. The wavelength and velocities are in the rest frame of the galaxy. 



Wavelength 

(A) 


FWHM 

(A) 


Norm 
(lO^'phs- 1 ) 


RMS Velocity 
(kins' 1 ) 


Flow Velocity 
(kms- 1 ) 


Transition 


21.6712;" 

qq 79+O.29 
OO. 1 ^_ .3i 


0.63l°; 3 2 2 
0.80l°;f 8 


11 0+ 21 

11 - u -2.0 

4 0+ 2 - 4 

^• u -2.3 


8700™ 
71001!°™ 


_i 40 +1800 


OVII triplet 
CVI Lya 


18.821°;!° 

99 1 9+0.13 
ZZ.1Z_ 12 


n +0.33 
U -0 

1H+ ' 26 

u - i,3 -0.13 


i -,+0.5 
L - 1 --0.4 

1 4+°- 7 


n +3200 
U -0 

1 7no +250 ° 

1 mu_ 1900 


-2400l 2 6 5 s °, 
+360lf 1 ° 


OVIII Lya 
OVII(f) 



3.3. Narrow Emission Lines 

Fig. [21 indicates that there are a number of significant 
narrow emission lines present in the spectrum. The two 
emission features either side of the CVI absorption line 
have been identified as a broad CVI emission line, and 
the 21. 7A, OVII(r), narrow emission line is confused by 
the presence of a broad O VII emission feature, so can no 
longer be reliably characterised. The results of a fit with 
two broad (O VII and C VI) and two narrow (O VIII and 
O VII(f)) emission lines are shown in Tables [21 and 

The fitted emission lines have shifted wavelengths rela- 
tive to the systemic velocity. These can be compared with 
those calculated from the line-finding program, without 
including any broad line features (Table If our identifi- 
cation of the emission lines is correct, the fit implies a flow 
velocity of approximately +360±400kms~ 1 (Tabled for 
the OVII(f) narrow emission line, at variance with the 
directly measured flow of -170±480kms" 1 (TableU}. The 
O VIII narrow emission line is found to have a flow veloc- 
ity of -2400±500kms~ 1 , in agreement with the directly 
measured results. 

3.4. Absorption Profile 

There are a number of spectral components available in 
SPEX that can now be added to the model to investi- 
gate the absorption lines that are present. The first model 
used here to analyse absorption is SLAB, which takes in- 



dividual ionic column densities, an outflow velocity and 
a RMS velocity and computes the absorption profile for 
each ion. SPEX varies the parameters to search for the best 
fit to the data. The outputs are then the best fit values 
for the column density, outflow velocity and RMS veloc- 
ity for each ion. These velocities give an initial indication 
as to whether the entire warm absorber is kinematically 
associated or whether there are separate regions moving 
with different velocities. 

In addition to the power law, two narrow and two 
broad Gaussian emission lines, we fitted individual ab- 
sorption profiles for CVI, N VII, OVII, OVIII, NelX and 
NeX (Table EJl, to our data. These ions were included be- 
cause they pro duce the most sig nificant lines in the spec- 
trum (Fig. |2J). iKriss et all l|2000|) find that one phase in 
their UV absorber has an ionization parameter log £= 1.06, 
and column density 5xl0 20 cm~ 2 , comparable with an X- 
ray absorbing gas. Therefore, we used the velocity profile 
of this phase as a starting point for our analysis setting the 
RMS velocity to that found for this UV phase: 40kms _1 . 

Some of the flow velocities are not well constrained 
making a detailed analysis difficult. However, we have 
plotted the velocities against the ionization parameter at 
maximum abundance to see if there is any connection be- 
tween the two and to see how they relate to the UV re- 
sults (Fig. |5J). From the slab results, in Table 01 we see 
that there appear to be two separate regions with different 
flow velocities, one at approximately -230 km s -1 and one 
at the systemic velocity. 
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Table 4. The parameters of the individually fitted SLAB 
absorption components, using a fixed RMS velocity value 
of 40kms~ 1 , with la error bounds. The flow velocity is in 
the rest frame of the galaxy. 



Species 


Column Density 


Flow Velocity 




logN (cm -2 ) 


(kms- 1 ) 


CVI 


17.38±g:i| 


-230+ 150 


NVII 


17.98lS:g? 


zw -180 


OVII 


17.801°;% 


-210+ 110 


OVIII 


18.13±H ? 


iq+120 
lu -110 


NelX 


17.401°;}° 


n+650 
u -360 


NeX 


17.63t»; 2 2 3 7 


fin +no 

DU -820 



An alternative SPEX model we use is XABS. This pro- 
vides a more physical view of the gas by providing a self- 
consistent model that contains all of the ions, in contrast 
to fitting each ion separately with SLAB. XABS applies ab- 
sorption by photoionized gas, with variable column den- 
sity, ionization parameter, elemental abundances, RMS 
velocity and flow velocity, to the model. The ionization 



parameter determines which transitions are likely to oc- 
cur in a given phase. By using this model we are able to 
see if there are multiple phases in the gas and ascertain 
the parameters for each phase. 

The data were fitted with different numbers of XABS 
models, the x 2 values for which are shown in Table 
The abundances of all elements were set to the solar val- 
ues. We see from the table the improvement in the quality 
of the fit found by adding multiple XABS components to 
the model. We conclude that three XABS components pro- 
vide a good match to the data. The best fit parameters 
(TableEJ) suggest three distinct regions, each with a differ- 
ent ionization parameter and outflow velocity. By plotting 
each phase of the model separately (Fig. 0) we can say 
which phase makes the most significant contribution to 
each absorption line. From the best fit parameters we see 
that there is a cool, high velocity phase producing the low 
ionization, log£=0.89, iron transitions forming the UTA 
at 16A, along with the NVI and OVI absorption lines. 
The intermediate phase, log£=2.14, produces the CVI, 
NVII, NelX and OVII absorption lines, whilst the low 
velocity, high ionization, log£;=3.26, phase produces the 
NeX, OVIII and highly ionized iron absorption features. 
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| 16 16.5 17 17.5 18 18.5 19 19.5 20 20.5 




31 31.5 32 32.5 33 33.5 34 34.5 35 35.5 36 

Wavelength (A) 

Fig. 4. A large scale view of Figure 01 the fit includes a power law, three XABS components, two broad Gaussians for 
the broad O VII and C VI emission lines and two narrow Gaussians for the narrow O VIII and O VII(f) emission lines. 

Table 5. The XABS best fit parameters for a model including a power law, two broad and two narrow Gaussian 
emission lines and three XABS components, with the 2a error bounds. The velocities are in the rest frame of the 
galaxy. 



Phase 


Flow Velocity 
(krns" 1 ) 


RMS Velocity 
(kmr 1 ) 


(erg cms -1 ) 


Column Density 
(10 21 cm- 2 ) 


Prominent Ions 


Phase 1 




n +20 


0.8912:3 


0-7912:11 


OVI, NVI, FelX-FeXIII 


Phase 2 


-260^0 


70+ 70 

' U -30 


2.141S;" 


0.7512:1? 


OVII, NVII, CVI, NelX 


Phase 3 


-601™ 


+30 


3.26l2i 8 7 


5.5lj 4 


OVIII, NeX, FeXVII-FeXX 
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SLAB model results 

□ XABS model results 

x Kraerner et al. 

a Kriss et al. 



onization (erg cm s 

Fig. 5. The relationship between the ionization parame- 
ter at maximum abundance and flow velocity of the el- 
ements in our SLAB and XABS models . Shown also are 
the UV results from iKriss et alJ l|2000|) (with conserva- 
tive tmcerta^ntie^^timated from their spectral resolution) 
and lKraemer et ahl l|2003f) including their given errors. 




Wavelength (A) 



Fig. 6. A model plot showing the contribution to the spec- 
trum by each of the three phases, the parameters of which 
are shown in Table Phase 1 is shown in black, Phase 2 
in dark grey and Phase 3 in light grey. 



Figs. El and 01 show the best fit, which includes a power 
law, two broad emission lines, at 21. 8A and 33. 7 A, two 
narrow emission lines, at 18. 9A and 22. lA, and three XABS 
components with the parameters listed in Table In ad- 
dition to the elements mentioned above there may also 
be contributions to the absorption profile from Silicon, 
Magnesium, and Argon. However, their impact on the 
spectrum is so small that they cannot be discussed in de- 
tail in this analysis. 

3.5. Ionization Structure 

Using the best fit values obtained in Section l3~H we are able 
to determine whether the three phases of the absorber are 
in pressure equilibrium (Krolik et al. 1981). In order to do 




Fig. 7. Thermal stability showing the equilibrium tem- 
p erature as a functi on of 5 for the SED described 
in lYaaoob et aT]l)2003|) . The positions relating to the three 
warm absorber phases are indicated along with their er- 
rors, displayed as the thick lines. For ease of compari- 
son the errors have also been indicated with vertical lines; 
Phase 1 shown with black dashed lines, Phase 2 with dark 
grey dotted lines and Phase 3 with light grey solid lines. 

this we use the pressure form of the ionization parameter 
(S), which is related to £ by 



L 



A-kct 2 P 



0.961 x 10 4 £ 
T 



(2) 



where L is the luminosity, r is the distance of the absorber 
from the ionizing source, P is the pressure and T the tem- 
perature (Steenbrugge et al. 2005). 

There is a specific range in ionization for which pres- 
sure equilibrium is possible (Krolik et al. 2001). In order 
to investigate this we used the X S TAR i onization code and 
the SED shown in lYaaoob et all l)2003|) to produce Fig. EJ 
Regions where dT/dB > are stable against isobaric ther- 
mal perturbations, with regions where dT/d'E. < 0, i.e. 
where the curve doubles back, being unstable. In order 
for the different phases to be in pressure equilibrium they 
must have the same value of S. Included in Fig [7] are 
the locations of the three phases of the warm absorber 
in Markarian 509. From this it is possible to see that the 
errors in 3 on Phase 1 and 3 marginally overlap. However, 
if Phase 3 had a compatible S value with Phase 1 then it 
would be on the unstable branch and so would be unlikely 
to remain there. Physically Phases 1 and 2 are more likely 
to be in pressure equilibrium because they are more sta- 
ble and have similar values of S. Therefore, it is unlikely 
that all three phases identified here coexist as part of the 
same structure, but Phases 1 and 2 could be in pressure 
equilibrium. 

4. Discussion 

4.1. Broad Emission Lines 

We have identified, and included in the fits, two broad 
emission lines due to O VII at 21. 7A, and C VI at 33. 7A, in 
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the RGS spectrum of Markarian 509. By fitting the lines 
with SPEX an RMS velocity value of 8700±2000kms" 1 
was found for the O VII line and 7100±4000kms~ 1 for the 
C VI line. We also identify a possible C VI RRC at 25. 2A, 
similar to that seen in NGC4051 (Ogle et al. 2004), com- 
patible with the CVI broad line. There is also a broad 
peak at 19. 9A: possible explanations for this include a 
N VII absorption edge (18. 6 A) or an Fe XXIII broad emis- 
sion line (19. 9A). However, for the NVII absorption edge 
to be large enough to account for the rise in the continuum 
the two fitted N VII absorption lines would have to be ap- 
proximately 50% deeper than is observed. Fe XXIII is also 
unlikely since we would expect to see other, stronger lines 
emitted from this species at 8.3A and lOA. The 25.2A and 
19. 9A features have not been included in the fits because 
they are not firmly identified with any strong expected 
transitions. 

The RMS velocity values of the C VI broad emission 
line suggests that it does not originate in any phase of the 
warm absorber, but is likely to arise in the broad line re- 
gion of the AGN close to the centre of the system. We have 
not given a flow velocity value for the OVII broad emis- 
sion line as it is likely to be a combination of the three lines 
of the O VII triplet making any velocity measurement un- 
reliable. This will also alter the fitted RMS velocity value; 
we only fit one broad Gaussian to the triplet. Therefore, 
each of the three emission lines is likely to have a smaller 
RMS value than that shown, but must still be broad since 
we only see one line. 

Broad emission lines, similar to those found here, have 
been seen in the soft X-ray spectra of other Seyfert 1 galax- 
ies such as NGC 5548 (Steenbrugge et al. 2005), Markarian 
279 (Costantini et al. 2005) and 3C273 (Page et al. in 
prep.). Steenbrugge et al. identify both the OVII, 21. 7A, 
and CVI, 33. 7A, broad emission lines, which are also our 
two most reliable broad lines. The lines in NGC 5548 
are fitted with a FWHM velocity value of 8 000 km s" 1 
taken from the corresponding UV broad lines. iKriss et alJ 
(200(J report broad OVI and CIII emission lines with a 
FWHM of 11000 kms" 1 in the UV spectrum of Markarian 
509, a good match with our broad X-ray emission lines. 

4.2. Narrow Emission Lines 

If we identify the 18. 8A line as an OVIII emission line, 
then it has an outflow velocity of 2400±500kms~ 1 . This is 
a much larger flow velocity than any other narrow emis- 
sion line in the X-ray or UV spectra of Markarian 509. 
However, the OVIII absorption line, at 18. 9A, is likely 
to remove the long wavelength side of the emission line 
making it narrower and moving the peak to an apparently 
shorter wavelength. We find this a more probable answer 
than identifying the line as an inflowing, high order, OVII 
transition with a flow velocity of +3000 kms -1 . 

The second narrow emission line included in our fit is 
the forbidden O VII line at 22. lA. This is part of the OVII 
triplet and so other lines are to be expected, but these have 



been overshadowed by the broad O VII emission line cen- 
tred at 2 1.8 A. Note the apparent flow velocities for our two 
narrow emission lines are significantly different (Table |3j); 
this is likely to be due to both lines being confused by 
the presence of other emission and absorption lines and 
so it is still possible that these lines originate in the same 
region of the AGN and have the same velocity structure. 
This broad line contamination is also the likely cause of 
the discrepancy between the fitted (Table EJ and the di- 
rectly measured (Table flow velocities for the OVII(f) 
narrow emission line. 

4.3. Characteristics of the Absorber 

In the spectrum of Markarian 509 we find three sepa- 
rate absorption phases: Phase 1 contributes the Fe M-shcll 
UTA and the OVI and NVI absorption features, Phase 
2 contains absorption from OVII, NVII, CVI and NelX, 
and Phase 3 produces the OVIII, NeX and highly ionized 
iron absorption features (Table EJ- We also identify minor 
contributions to the spectrum from Silicon, Magnesium 
and Argon. The confidence intervals for our ionization pa- 
rameters have no overlap in log£, so we are looking at 
at least three distinct phases in the gas; however, a con- 
tinuous ionization distribution throughout the gas cannot 
be excluded. We have shown that it would be possible for 
Phases 1 and 2 to coexist in the same structure as they 
have compatible velocity profiles and may have the same 
value of the pressure form of the ionization parameter, so 
could be in pressure equilibrium. 

We have improved the reliability of our results by us- 
ing a number of different techniques to analyse the ab- 
sorption profile. Fig. [5] shows our SLAB and XABS results 
along with UV information from previous analyses. We 
can see that all of the individually fitted SLAB absorp- 
tion lines have the same flow velocity, within the errors, 
as their main contributing phase, found from the XABS 
analysis. These results are also backed up by the flow ve- 
locities derived from the line-finding program (Table . 
Some of the values shown are not well constrained due to 
the complexity of the spectrum and the frequent confu- 
sion between spectral features. This is especially evident 
in the case of the Ne IX and Ne X absorption lines, both of 
which are confused by the presence of highly ionized iron 
features, thereby decreasing the reliability of the fitted 
values. There is also the additional problem of absorption 
lines being produced in more than one phase (Fig. EJ) caus- 
ing the apparent line width to increase and the apparent 
wavelength of the line to change. 

IKriss et all l|2000|) and lKraemer et all l|2003[) identify 
carbon and oxygen as the main contributors to the UV ab- 
sorption spectrum and also find some nitrogen and silicon, 
all of which we have identified in the RGS X-ray spectrum. 
Interestingly, in Fig. [3] we see that a number of the UV 
absorbers identified have a similar outflow velocity to our 
Phase 1 and are also within the ionization range that this 
phase covers. We are able to predict the UV spectrum that 
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our phases would produce and we find that Phase 1 con- 
tributes significantly to the carbon, oxygen and nitrogen 
features in the UV spectrum. We have used the ratios of 
the fluxes of these lines and the velocity profiles to find 
which phases are compatible and conclude that our Phase 
1 is compatible with the UV Phase 2 from the Kraemer et 
al. and Kriss et al. analyses. Our Phase 2 also contributes 
to the UV spectrum but we have been unable to reach a 
definite conclusion about which UV phase it corresponds 
to. 

Kriss et al. identify one of their regions, Phase 5, as 
being compatible with an X-ray producing gas, having a 
column density of 5xl0 20 cm~ 2 , ionization of log £=1.06, 
an outflow velocity of 5 km s _1 and a FWHM of 40 km s _1 . 
This ionization and column density are compatible with 
our X-ray Phase 1, but the velocity profiles do not corre- 
spond so we conclude that this particular UV phase is not 
compatible with any of our X-ray phases. 

We can also look back to the previous a nalyses of this 
warm absorber and compare our findings. IPounds et alJ 
l)200l|) observed a weak Fe M-shell UTA region, a Ne IX 
absorption line and the O VII narrow emission line triplet 
in the first XMM-Newton RGS spectrum of Markarian 
509 taken on 25 November 2000. These features are ob- 
served in our spectrum apart from an identification of all 
three lines of the OVII triplet. The analysis of Pounds 
et al. finds no evidence of carbon or nitrogen absorption 
features and also appears t o indicate an unusu al ratio be- 
tween the triplet lines (see iMewe et alJ l|2003h for triplet 
line ratios). As we have shown, two of the triplet emis- 
sion lines are not detected after the inclusion of the broad 
OVII emission line. Pounds et al. identify two separate 
phases producing the RGS spectrum but do not give any 
values for the ionization parameters or column densities 
for us to compare with our analysis. 

lYaaoob et alJ l|2003ft analysed Markarian 509 Chandra 
data and fitted two absorption edges from OVII and 
O VIII. These features are not observed clearly in our data, 
which may be due to our better definition of the UTA re- 
gion, mimicking the OVII absorption edge. The OVIII 
edge at 14. 2A is extremely weak and so is not a promi- 
nent feature in our RGS data. Yaqoob et al. also observed 
no clear evidence of an Fe M-shell UTA centred around 
16A. However, there is good agreement in the identifica- 
tion of Ne IX and Ne X lines, as well as O VIII absorption, 
which are some of the strongest lines in our RGS spec- 
trum as well. These lines were fitted by Yaqoob et al. 
with an ionization parameter of log£ = 1.76^ ; 14 and col- 
umn density 2.06+ ; 45 x 10 21 cm -2 , which lie between our 
Phase 1 and Phase 2 in ionization; however, they fitted 
only one phase of absorbing gas to the data and so the 
results are expected to differ somewhat. There is a good 
velocity agreement with our Phase 2, as Yaqoob et al. fit 
an outflow velocity of 200±170kms _1 and use a curve of 
growth analysis to derive a velocity width of 100 km s -1 . 
There will be differences between Yaqoob et al.'s and our 
results due to their data being taken with a higher spec- 



tral resolution and lower signal to noise ratio. Moreover, 
Markarian 509 has a variable luminosity (Jackson et al. 
2000, Carone et al. 1996) with Ij2-wkeV typically in the 
range (1.3-2.6) xlO^ergs" 1 (Weaver et al. 2001), which 
may alter the spectral profile over time. The 0.5-10 keV 
unabsorbed flux during the first observation with XMM- 
Newton was 2.6x 10^ 11 ergcm -2 s _1 , a factor of 2.8 less 
than the 0.5-10keV flux of 7.3x 10~ n ergcm~ 2 s" 1 dur- 
ing the Chandra observation. The second XMM-Newton 
observation occurred just 6 days after the Chandra obser- 
vation with a 0.5-10keV flux of 6.3xl0 _11 ergcm~ 2 s _1 . 

5. Conclusion 

We have used the high resolution X-ray spectra obtained 
with the Reflection Grating Spectrometer onboard XMM- 
Newton to study the characteristics of the Seyfert 1 galaxy 
Markarian 509 and present the most detailed analysis to 
date of the X-ray absorbing and emitting photoionizcd gas 
present in this system. 

The continuum has been fitted with a simple power 
law. The analysis suggests the presence of a warm absorber 
consisting of three phases: Phase 1 contains the high out- 
flow velocity, low ionization OVI and NVI absorption 
lines along with an Fe M-shell UTA; Phase 2 contributes 
the medium ionization OVII, NVII, CVI and NelX ab- 
sorption features, and Phase 3 produces the low outflow 
velocity, high ionization OVIII, NeX and highly ionized 
iron absorption features. However, we cannot rule out the 
possibility of a continuous ionization change throughout 
one gas system. It is possible that Phases 1 and 2 may co- 
exist in pressure equilibrium, but it is unlikely that Phase 
3 is also part of the same structure. We have identified 
Phase 1 as being compatible with a UV phase identified 
in previous studies. 

In addition to the absorption features identified here, 
there is also evidence of emission in the spectrum of 
Markarian 509. We identify two narrow emission lines due 
to OVII at 22. lA, and OVIII at 18.8A. We also iden- 
tify two broad emission lines, due to OVII at 21.7A, 
and CVI at 33. 7A, which have an RMS velocity of 
8000±3000 km s - 1 consistent with originating in the broad 
line region. 
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